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Abstract 

We address the following question. Take the Constrained Minimal Supersym- 
metric Standard Model (CMSSM) with the two CP violating SUSY phases 
different from zero, and neglect the bound coming from the electric dipole 
moment (EDM) of the neutron: is it possible to fully account for CP viola- 
tion in the kaon and B systems using only the SUSY contributions with a 
vanishing CKM phase? We show that the BR(S ^sj) constraint, though 
CP conserving, forces a negative answer to the above question. This implies 
that even in the regions of the CMSSM where a cancellation of different con- 
tributions to the EDM allows for large SUSY phases, it is not possible to 
exploit the SUSY phases to fully account for observable CP violation. Hence 
to have sizeable SUSY contributions to CP violation, one needs new flavor 
structures in the sfermion mass matrices beyond the usual CKM matrix. 

13.25.Es, 13.25.Hw, 11.30.Er, 12.60.Jv 
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I. INTRODUCTION 



Since the initial work of Kobayashi and Maskawa, the Standard Model (SM) of elec- 
troweak interactions is known to be able to accommodate the experimentally observed CP- 
violation through a unique phase, Sckm, in the Cabibbo-Kobayashi-Maskawa mixing matrix 
(CKM). However, the available experimental information, namely Ek and e'/e, is not enough 
to establish this phase as the only source of CP-violation. 

Most of the extensions of the SM include new observable phases that may significantly 
modify the pattern of CP violation. Supersymmetry is, without a doubt, one of the most 
popular extensions of the SM. Indeed, in the minimal supersymmetric extension of the SM 
(MSSM), there are additional phases which can cause deviations from the predictions of the 
SM. After all possible rephasings of the parameters and fields, there remain at least two new 
physical phases in the MSSM Lagrangian. These phases can be chosen to be the phases of 
the Higgsino Dirac mass parameter (0^ = Arg[/i]) and the trilinear sfermion coupling to the 
Higgs, {(pAo = Arg[74o]) In fact, in the so-called Constrained Minimal Supersymmetric 
Standard Model (CMSSM), with strict universality at the Grand Unification scale, these are 
the only new phases present. 

It was soon realized that for most of the CMSSM parameter space, the experimental 
bounds on the electric dipole moments of the electron and neutron constrained (pAo,fj. to be 
at most O{10^^). Consequently these new supersymmetric phases have been taken to vanish 
exactly in most studies of CMSSM. 

However, in the last few years, the possibility of having non-zero SUSY phases has again 
attracted a great deal of attention. Several new mechanisms have been proposed to suppress 
EDMs below the experimental bounds while allowing SUSY phases 0{1). Methods of sup- 
pressing the EDMs consist of cancellation of various SUSY contributions among themselves 
0, non- universality of the soft breaking parameters at the unification scale and approxi- 
mately degenerate heavy sfermions for the first two generations @]. In a recent work we 
showed that, in a model with heavy sfermions of the first two generations and in the large 
tan/5 regime, Ek and Eb could receive very sizeable contributions from these new SUSY 
phases. Similar studies including a larger set of experimental constraints have reported 
the impossibility of such large supersymmetric contributions Q. 

In this work, we are going to complete our previous analysis with the inclusion of all 
the relevant constraints in a CMSSM scenario. In doing so we adopt a different perspec- 
tive. We will assume from the very beginning that both supersymmetric phases are 0{1), 
ignoring for the moment EDM bounds 0. In these conditions, and taking into account other 



In this paper we restrict our discussions to the CMSSM. If one relaxes some of the constraints 
of this model, for instance by allowing for large gluino mediated CP violation with non-universal 
soft SUSY breaking terms, then it might still be possible to have fully supersymmetric e and e'/e 

i- 

^ EDM cancellations may be obtained through non-trivial relative phases in the gaugino mass 
parameters (see for instance the third paper in Ref. Q). However, for the discussion of the present 
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CP-conserving constraints, we will analyze the effects on the low energy CP-violation ob- 
servables, especially Sk and Sb- It should be noticed that the model used in ||^ can be easily 
obtained as a limit of the CMSSM by decoupling the first two generations of squarks and 
neglecting the intergeneration mixing in the sfermion mass matrices. Hence, our results in 
the more general CMSSM will include this model as a limiting case. 

In the next section we study the new sources of flavor mixing present at the electroweak 
scale in any supersymmetric model. In section 3 we are going to analyze neutral meson 
mixing, i.e. K^-K^ and B^-B^ mixings, with large supersymmetric phases. Section 4 
will be devoted to the study of the branching ratio of the decay b —>■ sj. In section 5 
we will show the impact of the measured b —>■ s'-f branching ratio on the supersymmetric 
contributions to Ek and Eb- Section 6 will present our conclusions and in Appendices A 
and B we collect, respectively, the formulas for the integration of relevant RGE's and the 
different loop functions appearing in the text. 

II. FLAVOR MIXING IN THE CMSSM 

The issue of flavor changing neutral current (FCNC) interactions in the CMSSM has 
been widely investigated in the literature. For the completeness of the discussion, we briefly 
recall those properties which will be relevant for our analysis. 

The CMSSM is the simplest Supersymmetric structure we can build from the SM particle 
content. This model is completely defined once we specify the soft-supersymmetry breaking 
terms. These are taken to be strictly universal at some ultra-high energy scale, which we 
take to be the Grand Unification scale [Mgut)- 

i^Q)ij = i^u)ij = i^Dij = = ijn\)i3 = ^0 ^ij^ 

ml, = = ml (1) 
rrig = m^y = = mi/2. 

That is, there is a common mass for all the scalars, mg, a single gaugino mass, mi/2, and 
all the trilinear soft-breaking terms are directly proportional to the corresponding Yukawa 
couplings in the superpotential with a proportionality constant Aqc^'^^. 

Now, with the use of the Renormalization Group Equations (RGE) of the MSSM, as 
explained in Appendix A, we can obtain the whole Supersymmetric spectrum at the elec- 
troweak scale. All the Supersymmetric masses and mixings are then a function of ttiq, mi/2, 
Ao, (f)A, (l>fj. and tan/5. We require radiative symmetry breaking to fix \fi\ and \Bfi\ with 
tree-level Higgs potential 0. 



paper, no explicit mechanism for such a cancellation is needed 

^The RGE's of the MSSM have received a vast amount of attention in the literature. However, 
in most of the previous analyses the SUSY phases (pA and (/)^ are switched off. For this reason we 
prefer to give the relevant RGE's with non-vanishing SUSY phases in Appendix A. 
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It is important to notice that, even in a model with universal soft breaking terms at some 
high energy scale as the CMSSM, some off-diagonality in the squark mass matrices appears 
at the electroweak scale. Working on the basis where the squarks are rotated parallel to 
the quarks, the so-called Super CKM basis (SCKM), the squark mass matrix is not flavor 
diagonal at Mw- This is due to the fact that at Mqut there exist two non-trivial flavor 
structures, namely the two Yukawa matrices for the up and down quarks, which are not 
simultaneously diagonalizable. This implies that through RGE evolution some flavor mixing 
leaks into the sfermion mass matrices. In a general Supersymmetric model, the presence of 
new flavor structures in the soft breaking terms would generate large flavor mixing in the 
sfermion mass matrices. However, in the CMSSM, the two Yukawa matrices are the only 
source of flavor change. As always in the SCKM basis, any off-diagonal entry in the sfermion 
mass matrices at Mw will be necessarily proportional to a product of Yukawa couplings. 
The RGB's for the soft breaking terms are sets of linear equations, and thus, to match the 
correct quirality of the coupling, Yukawa couplings or tri-linear soft terms must enter the 
RGE in pairs, as we can see in Eqs. (A.1-A.3) in Appendix A. 

In fact, in the up (down) squark mass matrix the up (down) Yukawas will also be di- 
agonalized and so will mainly contribute to diagonal entries while off-diagonal entries will 
be due to the down (up) Yukawa matrix. This means, for instance, that in this model the 
off-diagonality in the M^^^ ^ matrix will roughly be c • YuYu- With c a proportionality factor 
that typically is, 

as expected from the loop factor and the running from Mqut to Mw- Nevertheless, we have 
to keep in mind that this is simply a typical estimate and the flnal value of c can suffer a 
sizeable variation depending on many other factors not present in this simple estimate. 

On the other hand, this has clear implications on the tan/3 dependence of these off- 
diagonal entries of the sfermion mass matrices. In the basis where the down Yukawa matrix 
is diagonal, we can write the up and down Yukawas as, 

MMz) = . Vi^M Mu , Yn{Mz) = M, (3) 

V 2MvK sm p \/2Mw cos p 

with VcKM the Cabibbo-Kobayashi-Maskawa mixing matrix and M^^d the diagonalized mass 
matrices for the quarks. We can see in this equation that, for tan/3 > 1, the up Yukawa 
matrix will maintain similar values when going to large tan /3. Hence, the off-diagonal entries 
in the down squarks mass matrix will be roughly stable with tan /3. In the up squark mass 
matrix we have the opposite situation and the tan/3 dependence is very strong. In this 
case the off-diagonal entries depend on the down Yukawa matrix that grows linearly with 
tan /? for large tan /?. This means that we can expect the flavor change in the up and down 
squark mass matrix to be similar when tan/3 ~ mt/m;, ~ 40. While for tan/3 ~ 2 the flavor 
change in the up mass matrix will be approximately (tan /3 = 40)^/ (tan /3 = 2)^ = 400 times 
smaller (see Appendix A for details). These points also apply to the left-right sub-matrices 
where again flavor changing entries will be due to the opposite isospin Yukawa matrix. In 
fact, this left-right sfermion mixing only appears after electroweak symmetry breaking. The 
expression for these matrices in the SCKM basis is. 
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Mi"/ = I^^VckmAUMz) - HMz)\e'^' cot (3 M^j (4) 
Mif = ^AUMz) - |^(Mz)|e*^nan/?Mrf (5) 

Then, these left-right mixings will have an additional suppression proportional to the mass 
of the corresponding right-handed quark ( remember that AjjVi ~ AoMjj). This is always 
true for all the generation changing entries that are produced by the A matrices. However, 
in the down mass matrix, this suppression can be partially compensated by a large value 
of tan (3 in the diagonal terms proportional to fi. These are all well-known facts in the 
different studies of FCNC processes in the framework of the CMSSM P,p!0[l and imply that 



flavor mixing is still dominantly given by the usual CKM mixing matrix in W-boson, charged 
Higgs and chargino vertices. 

In this work, we are especially interested in CP violating observables. Then we must also 
consider the presence of observable phases in the sfermion mass matrices. In the following 
we will take the CKM matrix exactly real to isolate pure effects of the new supersymmetric 



phases |jTT[|. The sfermion mass matrices contain several physical phases that give rise to 
CP violation phenomena. In particular, before RGE evolution, these phases (0^, 0^) are 
confined to the left-right part of the sfermion mass matrix while both the left-left, tuq, and 



w- 



right-right, mfj j-,, matrices are real and diagonal. However this is not true anymore at M] 
(pA leaks into the off-diagonal elements of these hermitian matrices through RGE evolution. 
From the explicit RGE in the MSSM, Eq.( [A.l|) , it is clear that this phase only enters the 
(mQ)jj evolution through the combinations {Ai}A\j)ij or [ADA^^)ij. At Mqut these matrices 
have a common phase, and so the combination (^^^) is exactly real. So, to the extent that 
the A matrices keep a uniform phase during RGE evolution, no phase will leak into the niq 
matrices. However, we can easily see from Eqs.(A.2) and (A. 3) that this is not the case, and 
different elements of the A matrices are renormalized differently. In this equation, we can 
see that only the terms involving two Yukawa and one A matrix can produce a mismatch in 
the phases. Moreover, these terms will only be important when there are no small Yukawas 
involved. Then, we can expect a mismatch only on the off-diagonal elements involving the 
third generation. Keeping this in mind, the general form of the mg matrix at Mw in terms 
of the initial conditions is, 

ml{M^) = vt^ml + ^\^^Al + ryjf^m?/, + (r^^J^^e^*- + ^^^^^^ e-^^-)A,my, (6) 

where the coefficients t] are 3x3 matrices with real numerical entries. In this expression 
we can see that the presence of imaginary parts will be linked to the non-symmetric part 
of the viq^^ matrices. As is clear from the mass matrices in Appendix A (Eqs. (A.5-A.7) 
and (A.12-A.14)), these non-symmetric parts of roq are always more that three orders of 
magnitude smaller than the corresponding symmetric parts. This means that, in the SCKM 
basis, the imaginary parts of any mass insertion are present only in one part per 2 — 3 x 10^, 
and are always associated with (3,i) MI, Eqs. (A.8-A.11) and (A.15-A.18). A very similar 
situation was also found by Bertolini and Vissani, in the CMSSM with vanishing susy phases 
for the leakage of Sckm [0)0 • So, we conclude that in the processes we will consider, we 



can take both M*^")^^. M^'^^^j;^ as real to a very good approximation. 
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In the following we will analyze the new effects of this model on indirect CP violation 
in K and B systems. In doing so, we will use both the exact vertex mixing method and 
the Mass Insertion (MI) approximation |T^. Notice that the MI approximation is extremely 



good in the case of the CMSSM where all the off-diagonal entries are sufficiently small. The 
size of these off-diagonal entries directly gives in the MI approximation the amount of flavor 
changing induced by the sfermion mass matrices. A possible exception may arise in the 
stop and sbottom sectors that, in any case, could be diagonalized to ensure the validity of 



the MI approximation |T^. As we will see in the next sections, this is frequently useful to 



understand the exact results obtained in the vertex mixing method. 



III. INDIRECT CP VIOLATION IN THE CMSSM 

In the SM neutral meson mixing arises at one loop through the well-known W-hox. 
However, in the CMSSM, there are new contributions to AF = 2 processes coming from 
boxes mediated by supersymmetric particles. These are charged Higgs boxes (H^), chargino 
boxes (x^) and gluino-neutralino boxes {g, x^)- The amount of the indirect CP violation in 
the neutral meson Ai system is measured by the well-known Em parameter: 

where AMj^^ is the Ai-Ai mass splitting. Em depends on the matrix elements of the AF = 2 
Hamiltonian, H^ff^^, which can be decomposed as 



'^efr' = -^r^iyuyt,?{c,{^i)Q,{^i) + cmqm + cmqm) (s) 



where the relevant four-fermion operators are given by 

Qi = dlYql ■ dli.ql 

Q3 = dlq^R ■ 4qR (9) 

with q = s,b for the K and i?-systems respectively and a, j3 as color indices. In the CMSSM, 
these are the only three operators present in the limit of vanishing rrid- 

At this point, we are going to divide our discussion into two parts. We analyze separately 
the effective operator Qi that preserves chirality along the fermionic line, and the operators 
Q2 and Q3 that change chirality along the fermionic line. As we will see below the flavor 
mixing in the sfermion mass matrix and the experimental constraints on both kinds of 
operators are very different. 

A. Chirality conserving transitions 

In Eq.(^, Qi is the only operator present that does not involve a chirality change in the 
fermionic line. With respect to the associated sfermion, no chirality change in the sfermion 
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propagator will be needed, and so, the suppression associated with left-right sfermion mixing 
can be avoided. In general, Ci(/io) can be decomposed as follows 



Ci(/io) = C: 



(10) 



The usual SM ly-box, where all the couplings are purely left-handed can only contribute 
to this effective operator. However, with Sckm = 0, does not contain any complex phase 
and hence cannot contribute to the imaginary part in em- In any case it will always be, in the 
CMSSM, the dominant contribution to AM_\4. Similarly, the charged Higgs contribution, 
depends on the same combination of CKM elements with no other CP-violating phase 
0. So it will not contribute to our CP violating observable. 

Gluino and neutralino contributions to Cf'^ are specifically supersymmetric. They in- 
volve the superpartners of quarks and gauge bosons. Here, the source of flavor mixing is not 
directly the usual CKM matrix. It is the presence of off-diagonal elements in the sfermion 
mass matrices, as discussed in the previous section. From the point of view of CP violation, 
we will always need a complex Wilson coefficient. In the SCKM basis all gluino vertices are 
flavor diagonal and real. This means that in the MI approximation, we need a complex mass 
insertion in one of the sfermion lines. As explained in the previous section, these MI are 
proportional to Yukawa couplings and real up to 1 part in 2 x 10'^. The complete expressions 
for the gluino contributions to AF = 2 processes in the MI approximation can be found 
in |jT5[. The bounds obtained there for the real and imaginary parts of the mass insertions 
required to saturate AM^ and Ek are. 



where M is an average squark mass. 

In the CMSSM, as we can see in Appendix A, these mass insertions are much smaller. 
In particular, the fact that the bound on AMk ,the real part of the MI, is satisfied implies 
that the imaginary parts are at least two orders of magnitude below the required value 
to saturate Ek ■ Hence, no sizeable contributions to Ek from gluino boxes are possible. 
The situation in B^-B^ mixing is completely analogous: assuming that the minimum phase 
from the mixing observable in the B-factories is around 0.1 radians, we would need an 
imaginary contribution not more than one order of magnitude below the real one. With 
the arguments given above, this is clearly out of reach for gluino boxes in the CMSSM. 
Neutralino contributions are generally smaller than gluino due to smaller couplings with the 
same source of flavor mixing. In fact, although neutralino vertices in the SCKM basis also 
involve the complex neutralino mixings, any imaginary part on this operator will only be 
due to a complex mass insertion. This can be seen in the explicit expressions in where all 
neutralino mixings in this operator appear in pairs with its complex conjugate counterpart. 

Finally, the charginos also contribute to Ci{My/)^. In this case, flavor mixing comes 
explicitly from the CKM mixing matrix, although off-diagonality in the sfermion mass 
matrix introduces a small additional source of flavor mixing. 




(11) 



(5^2)LL|sin(20LL) <3-10-=^ 
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i,j=l k,l=l aja'Y V ^td ^tq) 

where Zk = M?^/M^^, Si = M'^JM^^, and VaqG^"'''^^ represent the couphng of chargino, i, 
and squark k to the left-handed down quark, q. Finally a is an intermediate up-quark 
index associated with the factorization of the CKM mixing matrix. The expression for this 
coupling is then, 

G'-"^(r?icj,.-^=^rj5,Q,.) (13) 

where Tm and Tuji are 6x3 matrices such that the 6x6 unitary matrix Tu = {TmTuji} 
diagonalizes the up-squark mass matrix, TuM^tIj = diag{M?^, ...,M|^). Cr is one of the 
matrices that diagonalize the chargino mass matrix through a bi-unitary transformation 
C^M-Cl = diagiM^±,M^±), with, 

M- = ( Mw cosp\ . . 

^ \Mw sin (3 \fi\e"t'^ J ^ ^ 

From these equations it is clear that G^"'^^"^ will in general be complex, as both 0^ and (pA 
are present in the different mixing matrices. The loop function Yi{a,b,c,d) is given in Eq. 
(B.l) of Appendix B. 

The main part of in Eq([l2D will be given by pure CKM flavor mixing, neglecting the 
additional flavor mixing in the squark mass matrix |iy,|T6|. In this case, a = a' and 7 = 7', 
we have. 



)^(M^) = E E E YM, zi, s,, s,)] (15) 

i,j=l k,l=l «7 y^td^tg) 



But, taking into account that Yi{a,b,c,d) is symmetric under the exchange of any pair of 
arguments we have, 

Q{a,k)iQia,k)j*Q{^,l)i*Q{^,l)j Y^(^Zk, Zi, Si, Sj) = (16) 



i (^Qia,k)iQ(a,k)j*Q{'y,l)i*Q('y,l)j _^ Q{a,k)i* Q{a,k)j Q{y,l)iQ{-f,l)j*'^ Yi{Zk, Zi 



1 "Sj; Sj) 



and so Ci'^^, is exactly real This is not exactly true in the CMSSM, where there is 
additional flavor change in the sfermion mass matrices. Here, some imaginary parts appear 
in the in Eq (|l2l). Being associated to the size of intergenerational sfermion mixings, 
these imaginary parts will be maximal for large tan /?. In Fig(|I]) we show in a scatter plot 
the size of imaginary and real parts of Ci in the K system for a fixed value of tan (3 = 40. 
The region of susy parameters explored in this and all of the following scatter plots is 
50 GeV < mQ,mi/2, Aq < 500 GeV and < (f)A,4>ti < 27r. With these initial conditions we 
impose that all squarks are heavier than 100 GeV with the exception of the stops that, as 
the charginos, are only required to be above 80 GeV. Furthermore we impose the constraint 
from the b ^ decay. Notice that, as we will see later, this is a conservative attitude, in the 
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sense that other constraints that we do not impose could only make our conclusions stronger. 
Under these conditions, we can see here that, in the CMSSM, this Wilson coefficient is 
always real up to a part in 10^. Fig(0) is the equivalent plot for the case of B^-B^ mixing. 
Here, imaginary parts are relatively larger but, in any case, out of reach for the foreseen 
B-factories. 

Taking this into account, from the point of view of experimental interest, we will always 
neglect imaginary parts in the Wilson Coefficient Ci within the CMSSM. Notice that this 



would not apply in a general model with non- universality at the GUT scale ||T5|, and each 
particular model should be considered separately. 



B. Chirality changing transitions 

From the point of view of flavor change and CP violation, operators Q2 and Q3 are 
different from Qi. These two operators always involve a change in the chirality of the 
external quarks, and consequently also a change of the chirality of the associated squarks or 
gauginos. In particular, this implies the direct involvement of the supersymmetric phases. 
On the other hand, these operators are suppressed by the presence of down quark Yukawa 
couplings, and so can only be relevant in the region of large tan/5 We can write the 
different contributions to C2 and C3 as, 

C2{Mw) = C^{Mw) + Cl{Mw) (17) 
C3(Mh/) = Ct\Mw) + C^{Mw) 

In flrst place, charged Higgs contributes only to C2, but parallel to the discussion for CY'^ ■, 
the absence of phases prevents it from contributing to Em- 

Gluino and neutralino boxes contribute both to Q2 and Qs- However flavor change will 
be given in this case by an off-diagonal left-right mass insertion. In the CMSSM these 
MI are always proportional to the mass and are never enhanced by large tan/3 values (see 
Eq.(^) of the right handed squark. This implies that these left-right flavor transitions from 
gluino will always be smaller in the CMSSM than the corresponding chargino contributions, 
where flavor change is directly given by the CKM matrix. In fact, this is already well-known 



for the case of 6 ^ 57 decay [|T7l, which is completely equivalent from the point of view of 
flavor change. 

Hence, the most important contribution, especially for light stop and chargino, will be 
the chargino box. Before the inclusion of QCD effects, it contributes solely to the coefficient 

2 6 T/* T/ T/* T/ ™2 

CAMw)-!^ ^ ^ {V;,V,,r 2M^cos^/3'' G G H 

Y2{zk,zi,Si,Sj) (18) 



i,j=l k,l=l a-ya'y' 



where Tnq/{\/2Mw cos P) ■ Vag ■ H'^"'^^^ is the coupling of chargino, and squark. A;, to the 
right-handed down quark g, with. 
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and 1^2(0, b, c, d) is given in Eq. (B.2). Unlike the Ci Wilson coefficient, C3 is complex even 
in the absence of intergenerational mixing in the sfermion mass matrices |p. In fact, the 
presence of flavor violating entries in the up-squark mass matrix hardly modifies the results 



obtained in their absence |10| , p!6| . So, in these conditions we have, 

Ci{Mw)=Yl E [Fs{3,k,3,l,t,j)-2F,{3,k,l,l,z,j) + F,{l,l,l,l,t,j)] (20) 

i,j=l k,l=3,6 

2 

2M^ COS"' p 

where we have used CKM unitarity and degeneracy of the first two generations of squarks. 
Due to the differences between H and G couplings, this contribution is always complex in 
the presence of susy phases. The most relevant feature of Eqs. (|I^) and (^Dj) is the explicit 
presence of the external quark Yukawa coupling squared, m"^/ {2M^ cos"^ P). This is the 
reason why this contribution is usually neglected in the literature P,P|JTB|. However, as we 
showed in |0, this contribution could be relevant in the large tan/5 regime. For instance, in 
B^-B^ mixing we have m^/(2M^ cos^ /3) that for tan/3 > 25 is larger than 1 and so, it is 
not suppressed at all when compared with the Ci Wilson Coefficient. This means that this 
contribution can be very important in the large tan (3 regime and could have observable 
effects in CP violation experiments in the new B-factories. However, in our previous work 
0, we did not include the additional constraints coming from 6 ^ 57 decay. In the next 
sections we will analyze the relation of Em with this decay, and the constraints imposed by 
its experimental measure. 



IV. S ^ 5*7 IN THE CMSSM 

The decay 6 — > 57 has already been extensively studied in the context of the CMSSM 
with vanishing susy phases fl^ . Because the branching ratio is a CP conserving observable. 



the presence of new phases will not modify the main features found in concerning the 
relative importance of the different contributions. However, in the presence of the new susy 
phases, these contributions will have different phases and will be observable through the 
interference. As we will see next, the experimental constraints will also have a large impact 
on the imaginary parts of the decay amplitudes. 

This decay is described by the following AF = 1 effective Hamiltonian 



nffr = -^VtlVt, E (21) 



where the relevant operators are given by 



Q2 = SL7f,CLCL7^bL, (22) 
Qr = ^,s,a^-^F,^bn, (23) 

Qs = ^-^SLcr^'G.M. (24) 
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Here C2(/io) = 1, and the Wilson coefficients Cj^s can be decomposed according with the 
particles in the loop, 



C-jiMw) = CY{Mw) + Ct/{Mw) + Cf{Mw) + "(Mv^) (25) 
C^{Mw) = Cf{Mw) + Ci{Mw) + ct{Mw) + Cf^'iMw) 

Among these contributions, the W penguin is exactly the same as in the SM and it does 
not depend on any supersymmetric parameters, it is simply a function of SM couplings and 
masses. This contribution is 

CY{Mw) = -\xt[QuF^{xt) + F2{xt)) 

Cf{Mw) = -\xtF^{xt) (26) 

with Xt = m'l/M'^ and Qu the charge of the up quarks. Similarly, in the charged Higgs 
penguins all the variables are known with the exception of M^- Again this contribution is 
unchanged by the inclusion of the new susy phases, 

C^{Mw) = -^{coe f3{QuF,{xt/xh) + F^ixt/xh)) + QuF^{xt/xh) + Fi{xt/xh)) 
C^iMw) = -^{cot^l3F,{xt/xh) + F^{xt/xrS) (27) 

where Xh = M^/M^. This contribution gives a sizeable correction to the b ^ s'j decay 
that constrains the mass of the charged Higgs in two Higgs doublet models or in the MSSM 
with low tan/3. However, in the case of moderate-large tan/3, chargino contributions may 



partially compensate this charged Higgs contribution relaxing the constraints IT^ • 

In addition to the and charged Higgs contributions analyzed above, there are three 
specifically supersymmetric contributions mediated by gluino, neutralino and chargino. In 
gluino or neutralino diagrams flavor change is due to the off-diagonality in the sdown mass 
matrix. Being left-right flavor off-diagonal transitions, they are suppressed by the mass 
of the b quark. Indeed, smallness of gluino and neutralino contributions has already been 
established in [|1^ where it was shown that, in the CMSSM, such contributions are roughly 
one order of magnitude smaller than the chargino contribution. 

Together with the and charged Higgs, the most important supersymmetric contri- 
bution will be, especially in the large-moderate tan /3 regime, the chargino contribution. In 
the W and charged Higgs contributions, the necessary chirahty flip for the dipole amplitude 
is always proportional to mi,. However, in the chargino penguin the chirality flip can be 
made either through a chargino mass insertion in the loop or through an external leg mass 
insertion proportional to rrif,. In fact, as pointed out in ||^, this enhancement due to m^i/rrib 
is partially compensated by the presence of the b Yukawa coupling. Nevertheless, this com- 
pensation is only effective for low values of tan/3. In terms of the chargino-quark-squark 
couplings used in the previous section, these contributions are, 

k=l i=l a,l3=u,c,t ^thVls 
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\f2Mw cos /5 

fc = l j = l Q-,/3=M,C,t ^th^ts 



_Hi-^k)iQ*iPM)iM2<LFl{z,, s.)) (28) 



\/2Mw cos /5 

with the loop functions defined in Appendix B. Similarly to the situation for the Wilson Co- 
efficient C3, we can, to a very good approximation, neglect the presence of intergenerational 
mixing in the up-squark mass matrix |]9|,p!0[|, then, 

fc=3,6i=l 



m^z mi, 



mb y/2Mwcos(3 

ctiM^) = E E {G^'^'^'G*^'''^'Fl{z,,s,) - G^'^'^^G*^'''^'Fliz„s,) + 

k=3,6 i=l 

^ ^J' - ..))) (29) 

rUb \/ 2Mw cos p ^ 

where, once more, we use CKM unitarity and degeneracy of the first two generations of 
squarks. 

The second term in Cy^g in Eq. (|5(]|), which corresponds to the chargino mass insertion 
in the loop, is dominant in the large tan/? regime. Notice that both and 
are products of the squark and chargino mixing matrices that can be OiV) (in the case of 
fiavor-diagonal stop mixings). Then, for stop and chargino masses around the electroweak 
scale, this term has an extra enhancement of 1/ cos/3. This means that, for large tan/?, we 
can approximate these Wilson Coefficients as, 

ct^M^) = E E^ ^J' J i/(^'^)-G-^^-^^'F^(..,..) -iJ(^'^^'G*^^'^^-F^(.i,..)) 

k=3,6i=l V2Mw COS (3 

kt^^et^i mb V2MwcosP 



V. B ^ 57 AND EM •• CORRELATED ANALYSIS 

As we have seen in Section III, chargino contribution to the C3 Wilson coefficient, Eq.(pO|), 
is the main contribution to indirect CP violation of the new supersymmetric phases for large 
values of tan /3. However, if we compare this Wilson coefficient with the chargino contribution 
to the decay b — >• S7, Eqs. (|20| ) and (^), we can see that both chargino contributions are 
deeply related. In fact, if we make a rough approximation and assume that the two different 
loop functions involved are of the same order, i.e.. 
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Y2{zk, zi, Si, Sj) ^ Fl{zk, Si) Fl{zi, sj) (31) 

we would obtain, 

2 

C,{Mw) = {C,{Mw)f-^ (32) 

Of course, this cannot be considered as a good approximation. As we can see from their 
exphcit expressions in Appendix B, the loop functions are clearly different. Anyway, they 
can be expected to give results of the same order of magnitude. So, the order of magnitude 
of C3 is determined by the allowed values of C7, as we will explicitly show below. 



To reach this goal, we will follow [|18|, where they constrain in a model-independent way 
new physics contributions to the Wilson coefficients involved in the 6 — S7 decay. In terms 
of these Wilson coefficients, the branching ratio, BR(S — is. 



BR(5 Xs-i) ~ 1.258 + 0.382|^7|2 + 0.015|^8p + 
1.395i?e[^7] + 0.161i?e[^8] + 0.083i?e[^7^8* 



+ (33) 



where = Ca{Mw) / {Mw). The different coefficients appearing in Eg . (^3]) are the SM 
renormalization group evolved contributions, that must be recovered in the limit = 1. The 



numerical values are taken from |^. We have not taken into account the errors associated 
with the choice of the scale and the restrictions on the photon energy that do not modify our 
conclusions. Now, using the experimental measure, BR(i? ^^7) = (3.14 ± 0.48) x 10~^, 
we can constrain the allowed values of the complex variables ^7 and ^g. In fact, we can 
already see from Eq. ([33|), that in the approximation ^7 ^ this is simply the equation of 
an ellipse in the Re[^-j\-Im[S,-j\ plane. In the case of supersymmetry with large tan/?, the 
new physics contribution to ^7 and ^g will be mainly due to the chargino. The allowed values 
of ^7 directly constrain then the chargino contributions to Cj{Mw) and indirectly constrain 
the values of C^{Mw)- 

In figure H, we show a scatter plot of the allowed values of Re^Cj) versus Im^C^) in the 
CMSSM for a fixed value of tan/5 with the constraints from Eq. (|35D . Notice, that a relatively 
large value of tan/?, for example tan/? > 10, is needed to compensate the W and charged 
Higgs contributions and cover the whole allowed area with positive and negative values. 
However, the shape of the plot is clearly independent of tan /?, only the number of allowed 
points and its location in the allowed area depend on the value considered. In this figure 
we take tan /3 = 40 because only a large value could give rise to observable CP violation . 
The values of C7 and Cg used here are the values obtained in the CMSSM for a given set of 
initial conditions. Although we do not use the approximation ,^7 ^ ^g this does not modify 
the elliptic shape of the plot. 

Figure ^ shows the allowed values for a re-scaled Wilson coefficient 6*3 (Mvk) = 
M^/mgCs{Mw) corresponding to the same allowed points of the susy parameter space 
in figure|^ . As we anticipated previously, the allowed values for C3 are close to the square 
of the values of C7 in figure |^ slightly scaled by different values of the loop functions. This 
is the proof of the importance of the 6 — 57 constraint on the chargino contributions to 
indirect CP violation. 
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We can immediately translate this result to a constraint on the size of the chargino 
contributions to e^. 



In this expression M_m, AM_m and F_m denote the mass, mass difference and decay constant 
of the neutral meson Ai^. The coefficient rj^lfi) = 2.93 [|1^ includes the RGE effects from 



Mw to the meson mass scale, /i, and B^{ij) is the B-parameter associated with the matrix 



element of the Q3 operator |T9 



Then, for the K system, using the experimentally measured value of AMk we obtain, 

2 

el = 1.7 X W-^—^Im[Cs] ^ 0.4 x lO^^/mfCs] (35) 

Given the allowed values of in Fig.^, this means that in the CMSSM, even with large 
susy phases, chargino cannot produce a sizeable contribution to Ek- We have seen in section 
III that also gluino and neutralino give negligible contributions in the CMSSM or in a model 
without off-diagonal soft-breaking terms at the GUT scale. Hence indirect CP violation 
in the kaon system will be mainly given by the usual SM box and the presence of a CP 
violating phase in the CKM matrix, 5ckm is still needed. 

The case of B^-B^ mixing has a particular interest due to the arrival of new data from 
the B-factories. In fact, as explained at the end of section III and in 0, in the large tan/3 
regime chargino contributions to indirect CP violation can be very important. However, for 
any value of tan (3 we must satisfy the bounds from the b ^ s'y decay. So, if we apply these 
constraints to the B^-B^ mixing, 

2 

4 = o.l7^/m[C3] ^ 0.5 X 10-^/m[C'3] (36) 

where once again, with the allowed values of Fig.^ we get a very small contribution to CP 
violation in the mixing. We must take into account that the mixing-induced CP phase, 
9m, measurable in B'^ CP asymmetries, is related to Eb 9m = arcsin {2V2 -eb}. The 
expected sensitivities on the CP phases at the B factories are around ±0.1 radians, so 
this super symmetric chargino contribution will be completely out of reach. Gluino and 
neutralino contributions to indirect CP violation can also be discarded in the CMSSM. 
Once again we have to conclude that no new contributions to indirect CP violation from the 
new susy phases will be observable in B^ CP asymmetries in the framework of the CMSSM. 
Recently, CDF has provided preliminary indications that sin 2/3 is in agreement with 



the SM predictions. Clearly, from the above result, Eq. (|3q) , it appears that the CMSSM 
contribution is too small by itself to account for this result. 



VI. CONCLUSIONS 

In this work, the effects of non- vanishing supersymmetric phases on indirect CP violation 
in K and B systems have been analyzed within the CMSSM. We have found that operators 
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involving only left-handed external quarks are not sensitive to these new phases at an 
observable level. This is due to the absence of intergenerational mixings beyond those 
originated from the CKM matrix. On the contrary, operators involving both right and left- 
handed quarks are in general complex, even in the absence of Sckm, and could be relevant in 
the large tan j3 regime. However, we have shown that these contributions are deeply related 
with the BK{B —>■ Xg'j) decay. So, taking into account the constraints coming from this 
decay these contributions also turn out to be too small to be measured experimentally. 

Although these conclusions are specific for indirect CP violation, they could also be 
implemented for chargino mediated direct CP violation in the decays. Again, in these 
decays the same chargino-quark-squark couplings are involved and we can also expect a 
big impact of the 6 — 57 constrain. In fact, the conclusions reached in this paper are 
far more general. The correlation between 6 — 57 and susy induced indirect CP violation 
exists in any supersymmetric model with sufficiently small intergenerational mixings in the 
sfermion mass matrices. This would include specifically all the models without new fiavor 
structures beyond the usual CKM matrix at the GUT scale and simplified models as the 
one the authors used in p. 

In summary, concerning the simpler supersymmetric models, like CMSSM, the con- 
straints coming from BR(i? — >■ X^j) decay are sufficient to rule out pure supersymmetric 
indirect CP violation in K and B systems, even in the absence of any electric dipole mo- 
ment constraints. This has very important consequences for the supergravity induced models 
where a cancellation between different supersymmetric contributions allows large supersym- 
metric phases while respecting EDM bounds 0. In these models, even in the regions of 
parameter space where this cancellation occurs, no observable effect of the large susy phases 
will appear on indirect CP violation experiments. However, as pointed out by Back and Ko 
0, these phases would still be observable in CP asymmetries in the 6 — 57 decay. 

All this means that the presence of large susy phases is not sufficient to produce ob- 
servable effects at the low energy experiments. In particular, new sources of flavor change 
beyond the usual CKM matrix are needed. And so, any deviation from the SM expectations 
at indirect CP violation experiments due to supersymmetry should be taken as a sign of 
non-universality of the soft breaking terms. In this context one recalls the recent studies on 



superstring compactifications with non-universal gaugino masses [21 
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APPENDIX A. INTEGRATION OF RGE'S IN CMSSM 



In this Appendix we describe the new features of the integration of RGE's in the CMSSM 
with non-vanishing susy phases relevant to our analysis. The complete matrix form of the 
RG equations can be found in 0. Using their notation and conventions, with the only 
change of Ag = niY^, we will mainly concentrate on the left-left scalar-quark mass matrix 
and the tri-linear soft breaking coupling evolution, 

dml _ 



dt 



YvY^ml + mlYuY^ + 2( Wl^^I + UlYuY^ + AuA^ + 



(A.l) 



I (fag + 3^2 + Au - (f 03^3 + SazMa + |aiMi) Yu - 

\[AAuY^Yu + QTr{AuY^)Yu + hYuY^Au + ?,Tr{YuY^)Au + 

2AdYIYu + YdYIAu] (A.2) 

I (fas + 3^2 + - (fasMg + SasMa + |aiMi) fz) - 

\[AAdYIYd + 6Tr(iz5Vi)l^D + ^YdY},Ad + 3rr(y},yi)iB + 

2i^Y-jF^ + YuY^Ad + 2Tr(isFi)r^ + Tr{YEYl)An\ (A.3) 

Except for the Yukawa coupling matrices, the RGE's of all other quantities are linear 0. 
This means, in particular, that RGE's of all soft masses, though coupled, can be solved as 
a linear combination of the GUT-scale parameters mo, A^e^^^ and M1/2 at any scale below 
Mg- However, one notices that the initial conditions on the trilinear couplings require 
the knowledge of the particular Yukawa texture at the unification scale. To do this, we 
numerically integrate the Yukawa RGE with a given value of tan/? and in terms of the 
fermion masses and the CKM matrix. Specifying the GUT-scale initial conditions in this 
way, it is straightforward to compute all soft masses at My/ for arbitrary values of mo, 
Aqc^'^^ and Mi/2- Thanks to the linearity of the corresponding RGE's the soft masses at 
Mw admit the following expansion, 

Au,d{Mz) = ^oe*^^ + «?/,dMi/2 
ml^uA^z) = ^Q,u,D ml + Vq,u,d + Vq,u,d ^1/2 (A.4) 

+ ylQ,U,D^ +Vq,U,D^ J ^0^^1/2 

where the coefficients a and 77 are 3x3 matrices with real numerical entries. One notices that 
the matrices rrtq jj ^{Mz) would be completely real were it not for the non-symmetric terms 
in the matrix f]Q^u However, it will be seen from the specific examples that this matrix 



dA 



u 



dt 



dA 



D 



dt 
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remains nearly symmetric, and thus, CP-violating entries mQjjjj{Mz) are extremely sup- 
pressed. Moreover, one notices that Au^dIMz) carries, in general, large CP violating phases; 
however, these terms are effective only for intragenerational Li?-type mixings. Hence, this 
particular observation shows the importance of chargino contributions for CP violation in 
FCNC processes, as explained in sec. II. 

As mentioned before, due to the non-linearity of the RGB's for Yukawa matrices, it is 
not possible to give a fully analytic solution for the soft mass parameters. Nevertheless, 
once we fix tan (3, we can numerically integrate the Yukawa RGE. Therefore, below we give 
semi-analytic solutions of RGE's for tan (3 = 2 and tan /5 = 40 to illustrate the small and 
large tan (3 regimes. 

Fixing tan/3 = 2, we get for the relevant rj matrices in Eq.( |A.4| ), 



q 



+ Vq ) 



(9 A) 



I 7.07 

2.79 X 10"^ 
\ -7.02 X 10-3 

5.34 X 10-6 
-3.44 X 10-5 
7.90 X 10-^ 






2.79 X 10"^ -7.02 X 10-^ 
7.07 4.92 X 10-2 
5.74 

7.90 X 10"^ 



4.92 X 10-2 

-3.44 X 10-5 
2.29 X 10"^ 
-5.52 X 10-3 





-5.52 X 10" 
0.15 

1.34 X 10-^ 







-1.34 X 10-s 8.55 X 10" 



-8.55 X 10- 




(A.5) 



(A.6) 



(A.7) 



where the vanishing off-diagonal entries in the last matrix mean values smaller than 10"^° 
in absolute magnitude. Among the matrices involved in Eq.( [A.4| ), r]^ is always the largest 
one for similar values of M1/2 and rriQ. So, it sets the scale of the matrix element while rj^sA) 
is the only one that can produce an imaginary part. Hence, we do not specify the other 
77-matrices, which are not important for our discussion. 

Once we obtain the mQ{Mw) matrix with the help of Eq.( |A.4| ) we can get the values of 

the and in the SCKM basis, that give the size of flavor change in the squark 

mass matrices compared with the diagonal elements. For tan/? = 2, those elements of the 
squark mass-squared matrix causing LL transitions between first and second, as well as 
second and third generations, are given by 



(M 



-2.79 X 10-^ m2 



9.30 X 10- 
+ 8.15x10-' AoMi/2 cos 0A 
= -4.07 X 10-5 m2 - 1.15 x 10" 



a2 

^0 



1.17 X 10-6 M- 



1/2 



/2 

+ 1 X lO-MoMi/2Cos0A - 1.71 X lO-MoMi/2 ^ sin 0a 



' Al - 1.61 X 10" 

-1-7 



Ml 



LL )l2 



9.38 X 10- 



rrin 



3.75 X 10-^2 + 2.79 x 10"^ 



/2 



+ 6.87 X 10-5 AoMi/2 COS0A 



23 



1.67 X 10 



-2 



rrin 



5.32 X 10"^ Al + 4.91 x IO-2 



1/2 



- 1.1 X 10-2 AoMi/2 cos 0A - 1-70 X lO-MoMi/2 isin0A 
Now, we repeat the same quantities above for tan/? = 40: 



(A.8) 
(A.9) 
(A.IO) 
(A.ll) 
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AgA) . (gA)T. 



/ 7.07 2.44 X 10"^ -5.80 x lO'^ 
2.44 X 10"^ 7.07 4.06 x lO'^ 

y -5.80 X 10-3 4.06 X 10-2 4.97 

.32 X 10-*^ -4.57 X 10-^ 



^.62 X iU " -4.5Y X iU " ( .^Z X 
-4.57 X 10-5 5.20 X 10-^ -5.47 
7.82 X 10-^ -5.47 x 10"=^ 0.: 



7.82 X 10^^ 
X 10-3 
0.22 

-1.64 X 10-6 \ 

1.14 X 10-5 

^ 1.64 X 10-6 -1.14 X 10-5 y 









(A.12) 
(A.13) 
(A.14) 



(MIl ) 12 = -8-77 X 10-5 - 2.77 x 10-5 - 3.0 x 10"^ 

+ 1.21 X lO-MoMi/2COS0A + il.l X lO-^MoMi/2 sin (A.15) 

(M|^^"^)23 = -1.28 X 10-2 - 2.70 x 10-^^ - 3.77 x 10"^ M^/^ 

+ 5.67 X lO-MoMi/2 cos 0A + ? 2.30 x 10"MoMi/2 sin 0^ (A. 16) 

{MIl'^^)i2 = 7.51 X 10^5 ^ 7 74 ^ IQ-^Al + 2.44 x M^/^ 

-9.13x10^5 AoMi/2 cos 0a (A. 17) 

= 1-34 X 10-2 + 7.84 x 10-^^^4.05 x 10^2^2^^ 

- 1.1 x 10-2 AoMi/2 cos (f)A + i 2.28 x 10-5 ^0^1/2 sin 0a . (A.18) 

A comparison of the corresponding quantities in tan P — 2 and tan P — AO cases reveals 
the sensitivity of the results on tan/3. As explained in Sec. II, Yu{Mz) remains nearly 
unchanged while Yr,[Mz) assumes an order of magnitude enhancement as tan/9 varies from 
2 to 40. This change in Yd{Mz) affects various quantities as dictated by the differential 
equations (A.l - A. 3). 
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APPENDIX B. LOOP FUNCTIONS 



In this appendix we collect the different loop function used in the text. The functions 
Yi and Y2 entering B-B and K-K mixings are given by 



Yi{a,b, c, d) 



Ina + 



{h — a){c — a){d — a) 
(a — c)(6 — c){d — c) 



a - b){c - b){d - b) 
d^ 



(a — d){b — d){c — d) 



In 6 
Ind 



(B.l) 



and 



Y2{a, b, c, d) 



Ucd 



+ 



{b — a)(c — a){d — a) 
c 



Ina + 



■ In c + 



{a-b){c-b){d-b) 
d 



In 6 



In d 



{a - c){b - c){d - c) {a - d){b - d){c - d) 

For the analysis of fo — > 57 branching ratio the following loop functions are relevant: 

1 



(B.2) 



F,{x) 
F2{x) 

FAx) 



i2{x - ly 
1 

i2(x - ly 
1 

2{x-lf 

1 

2{x - 1)3 
1 



{x-^ - + 3x + 2 + 6xlnx), 
{2x^ + 3a;^ - 6x + 1 - Inx), 
[x^ -4a; + 3 + 21na;), 
[x^ — 1 — 2xlnx), 



Fl{x, y) = - [QuF2{y/x) + F^{y/x)\ , 

X 

FU^,y) = ^ [QuF,{y/x) + Fs{y/x)] , 

Fl{x,y)^-F2{y/x), 

F^{x,y) = -F4{y/x). 

X 



(B.3 

(B.4 

(B.5 

(B.6 

(B.7 
(B.8 
(B.9 
(B.IO 
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FIG. 1. Imaginary and Real parts of the Wilson coefficient in Kaon mixing. 
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FIG. 2. Imaginary and Real parts of the Wilson coefficient in B mixing. 
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FIG. 3. Experimental constraints on the Wilson Coefficient C7 
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FIG. 4. Allowed values for the re-scaled WC C3 



